Astronomy & Astrophysics manuscript no. ms'final 


July 1, 2009 


(DOI: will be inserted by hand later) 





A multiwavelength study of Swift GRB 060111 B constraining the 
origin of its prompt optical emission * 

G. Stratta 1 , A. Pozanenko 2 , J-L. Atteia 3 , A. Klotz 4 , S. Basa 5 , B. Gendre 5 , F. Verrecchia 1 , M. Boer 6 , S. Cutini 1 , M. 
Henze 7 , S. Holland 8 , M. Ibrahimov 9 , F. Ienna 3 , 1. Khamitov 10 , S. Klose 7 , V. Rumyantsev 11 , V. Biryukov 12 , D. 

Sharapov 9 , F Vachier 13 , S. Arnouts 514 , and D.A. Perley 15 

1 ASDC, ASI Science Data Center, Via Galileo Galilei, 00044 Frascati, Italy 

2 Space Research Institude IKI, Profsoyznaya 84/32, 1 17910 Moscow, Russia 

3 LATT, Universite de Toulouse, CNRS,14 Av. E. Belin, F-31400 Toulouse, France 

4 CESR, Universite de Toulouse, CNRS, BP 4346, F-31028 - Toulouse Cedex 04, France 

5 LAM, Observatoire Astronomique de Marseille, Pole de l'Etoile Site de Chateau-Gombert 38, rue Frederic Joliot-Curie 
13388 Marseille cedex 13, France 

6 Observatoire de Haute Provence, F-04870 Saint Michel 1' Observatoire, France 

7 Thiieringer Landes stern warte Tautenburg, Sternwarte 5, D-07778 Tautenburg, Germany 

8 NASA, Goddard Space Flight Center, Greenbelt, MD 20771, USA 

9 Ulugh Beg Astronomical Institute, Tashkent 700052, Uzbekistan 

10 TUBITAK National Observatory, Akdeniz Universitesi, 07058, Antalya, Turkey 

11 SRI Crimean Astrophysical Observatory, Nauchny, Crimea, 98409, Ukraine 

12 Crimean Laboratory of the Sternberg Astronomical Institute, Nauchny, Crimea, 98409, Ukraine 

13 IMCCE, Observatoire de Paris, 77 Avenue Denfert Rochereau, F-75014 - Paris, France 

14 Canada France Hawaii Telescope corporation, 65-1238 Mamalahoa Hwy, Kamuela, HI 96743, USA 

15 Department of Astronomy, University of California, Berkeley, CA 94720-341 1, USA 

Received ; accepted 

Abstract. The detection of bright optical emission measured with good temporal resolution during the prompt phase of GRB 
0601 1 IB makes this GRB a rare event that is especially useful for constraining theories of the prompt emission. For this reason 
an extended multi-wavelength campaign was performed to further constrain the physical interpretation of the observations. In 
this work, we present the results obtained from our multi-wavelength campaign, as well as from the public Swift/BAT, XRT, 
and UVOT data. We identified the host galaxy at R ~ 25 mag from deep i?-band exposures taken 5 months after the trigger. Its 
featureless spectrum and brightness, as well as the non-detection of any associated supernova 16 days after the trigger, enabled 
us to constrain the distance scale of GRB 0601 1 IB within 0.4 < z < 3 in the most conservative case. The host galaxy spectral 
continuum is best fit with a redshift of z ~ 2, and other independent estimates converge to z ~ 1-2. From the analysis of the 
early afterglow SED, we find that non-negligible host galaxy dust extinction, in addition to the Galactic one, affects the observed 
flux in the optical regime. The extinction-corrected optical-to-gamma-ray SED during the prompt emission shows a flux density 
ratio F y jF opt = 10~ 2 - 10~ 4 with spectral index /3 r . op , > f3 y , strongly suggesting a separate origin of the optical and gamma-ray 
components. This result is supported by the lack of correlated behavior in the prompt emission light curves observed in the two 
energy domains. The temporal properties of the prompt optical emission observed during GRB 06011 IB and their similarities 
to other rapidly-observed events favor interpretation of this optical light as radiation from the reverse shock. Observations are in 
good agreement with theoretical expectations for a thick shell limit in slow cooling regime. The expected peak flux is consistent 
with the observed one corrected for the host extinction, likely indicating that the starting time of the TAROT observations is 
very near to or coincident with the peak time. The estimated fireball initial Lorentz factor is T > 260-360 at z = 1-2, similar to 
the Lorentz factors obtained from other GRBs. GRB 0601 1 IB is a rare case of a GRB with both a bright, well-observed optical 
counterpart and a 'canonical' early X-ray light curve, thus providing a good test case of the reverse shock emission mechanism 
in both energy ranges. 

Key words, gamma-ray burst - GRB 0601 1 IB - host galaxy - redshift 



Send offprint requests to: G. Stratta Palma jointly by Denmark, Finland, Iceland, Norway, and Sweden, in 

* Based on observations performed with: TAROT at the Calern ob- the Spanish Observatorio del Roque de los Muchachos of the Instituto 
servatory. The Nordic Optical Telescope, operated on the island of La de Astrofisica de Canarias. CFHT, which is operated by the National 



2 



G. Stratta et al.: A multiwavelength study of Swift GRB 0601 1 IB constraining the origin of its prompt optical emission 



1. Introduction 

Measuring the multi-wavelength spectrum of the prompt emis- 
sion of gamma-ray bursts (GRBs) is a challenge, requiring 
coordinated observations both from space and on the ground 
that are perfomed on very short notice (a few seconds). Such 
observations are nevertheless sorely needed to constrain the 
natur e of the prompt emission . Rapi d -response telescopes 



(e.g. Akerlof et al 



2004; Covinoetal 



2003 



2004 



Boerj 200ll Perez-Ramirez et al 



Bloo m et al.l |2006) and develop 



ments in the technique of wide-field optical surveys (e.g . 
Vestrand et al.l2002l;|Pozanenko et ap2004l; ICwiok et all2 005 ; 
Tamagawa et al.l I2005I) are fundamental to this purpose. So 
far, a few dozen GRBs have been detected at optical wave- 
leng ths when the GRB was still active, or soon after its en d 
(e.g. iBlake et alJliool lYuan et alJEooi iGendre et al.ll2008bh . 



In several cases (e.g. GRB 990123, GRB 0803 19B), the op- 
tical flux after the peak emission shows a very steep decay 
that renders the detection of this component even more chal- 
lenging (e.g. lMelandri et aD l2008). The demanding rapidity of 
observing these bright optical flashes (about a dozen seconds 
after the trigger or less) is on average fulfilled by automatic 
ground-based telescopes but not by the onboard UV- Telescope 
SwiftfWOT, whose first obse rvation times range from 40 to 
200 seconds after the trigger (Romin g et al.l 2009). It is thus 



difficult to infer the true fraction of GRBs showing this emis- 
sion. 

The common interpretation of the observed emission from 
GR Bs relie s on the fireball internal-external shock model (see 
e.g. IZhangl 2007 for a recent review). Whether the origin of 
the optical counterpart observed during the prompt emission is 
from internal shocks, like the high-energy emission, or from 
an external (revers e) shock, is still an open debate in the con- 
text of this model dGenet et al . 2007b) . Answering this question 
will provide useful insight into the formation of GRB shocks. 
In general, the early optical emission might be a superposition 
of internal shocks, reverse shock and forward shock, with dif- 
ferent contributions at different phases of the time history. 

The first Swift-BAT position notice of GRB 0601 1 IB was 
distributed at 20:16:03 UT, 20 s after the trigger, allowing 
the TAROT and ROTSE III robotic telescopes to observe 
during prompt em ission dKlotz et al.l 120061: lYost et aL 2006; 



Rvkoff et al.ll2009l) . TAROT Calern (France) started observing 
at 20: 16: 1 1 UT (28 s after trigger) and detected a bri ght optical 
emission with R~13.7, fading as r 2A±0A until 80s dKlotz et al .1 



2006h . The position of the afterglow as measured by MITSuME 
from 45-minute exposure images taken about 10 minutes af- 
ter the trigger is R.A.= 19* 05 ffl 42.4 7 \ Dec.= +70° 22' 33.1 " 
(J2000.0), with a precision of 0.3" dYanagisawa et al.ll2006l) . 



Research Council of Canada, the Institut National des Sciences de 
l'Univers of the CNRS of France, and the University of Hawaii. 
The 2.6m Shajn telescope at the Crimean Astrophysical Observatory. 
Some of the data presented here were obtained at the W.M. Keck 
Observatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Observatory 
was made possible by the generous financial support of the W.M. Keck 
Foundation. 



This is fully compatible with the position measured by UVOT 
with an accuracy of 0.5" (RA = 19 /, 05'"42.48 5 , Dec.= +70° 22' 
33.6"|PerrietaL|2006a). 

Because of the high declination of the source, ground-based 
optical follow-up has been challenging: the position cannot 
be observed from any major South Hemisphere observatory, 
and while it is circumpolar in the northern hemisphere, at late 
times the hour angle of the source is very large throughout 
the night and not accessible to equatorial telescope mounts. 
Despite these unfavorable conditions we set-up an observing 
campaign that provided many important results: the identifi- 
cation of the host galaxy, strong constraints on any supernova 
associated with the GRB and on the host-galaxy redshift. We 
present in Table|2]the list of the observations of the GRB error 
box that led to the results presented in this paper. 

The first part of this paper deals with the observations of 
GRB 0601 1 IB (§2). Results from the multiwavelength prompt 
and afterglow analysis are presented in $3] Results from the 
host galaxy identification, as well as the redshift constraints of 
the source, are presented in §|4] while in §0 we discuss the na- 
ture of the prompt optical emission in light of our multiwave- 
length campaign results. 

2. Observations and data analysis 

2.1. Swift/BAT and WAM/Suzaku 

GRB 06011 IB was a moderately bright GRB detected on 
2006, January 11 at 20:15:43.24 UT (her eafter W) by the 



BAT instrument o n the Swift spacecraft (Pet ri et al.l 12006a: 



Tuel ler et all [2006) 



This burst exhibits two peaks reaching 
their maxima at t tI j g +2 s and t tr i g +55 s. While the first peak 
is clearly detected in the full BAT energy band, the second 
is significantly softer (Fig. [TJ. The duration of the burst is 
= 59 ± 5 s in the range 15-150 keV. The time-averaged 
spectrum (from t tl -i g - 2 to t^g + 63 s) is well-fit by a power- 
law model, with best-fit photon index F r = 0.9 + 0.1. The ob- 
served peak energy flux and fluence in the 15-150 keV band are 
2.6 X 10~ 8 erg cm 2 s _1 and 1.6 x 10~ 6 erg cm~ 2 , respectively. 

The prompt high-energy emission was also observed by 
WAM on Suzaku, with a duration of Tgo = 25.4 s and a 
fluen ce of 5.6 ± 0.8 x 10~ 6 erg/cm 2 in the 100 - 1000 keV 
band (Sat o et al . 2006). A joint spectral analysis of the prompt 
emission with Swift/BAT and Suzaku/WAM data assuming a 
Band model, provide s a best-fit, observer-fr ame peak energy of 
Epeak = 462+ 2 ^ keV dFukazawa et al.ll2007l) . This value is high 
enough to consider the spectral index measured by Swift to be 
the low-energy spectral index of the Band function describing 
the energy spectrum of the prompt emission. 



2.2. Swift/XRT 

Swift s lewed immed iately to the burst, and the X-Ray Telescope 
(XRT iBurrows et al.l 2005) began data acquisition at t tr j g + 87 
s in Windowed Timing (WT) mode. At t tl i g + 153 s the count 
rate was at a critical level for XRT acquisition mode and an 
automatic flickering between the WT mode and the Photon 
counting (PC) mode started. At that time, the afterglow was too 
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Fig.l. Swift/BAT count rate of GRB 06011 IB in the ranges 
(from top to bottom): 15-25 keV, 25-50 keV, 50-100 keV, 100- 
150 keV and 15-150 keV. Count rate units are counts per sec- 
onds per illuminated-detector (note illuminated-detector = 0. 16 
cm 2 ). 
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Fig. 2. Swift/BAT count rate (in units of counts per seconds 
per illuminated-detector) of GRB 0601 1 IB in the range 15-150 
keV. The TAROT light curve superimposed on the BAT count 
rate is cr = £10~ 04m ' i with k = 10 4 in arbitrary units. The tem- 
poral ranges covered by UVOT and XRT data are also plotted. 



bright for the PC mode, and the PC data were affected by mod- 
erate pile-up. At t tr j g + 426 s, the XRT switched permanently to 
the PC mode and continued observations until t tr ; g +2.3 days. In 
the following, we used WT data taken between ttn g + 87 s and 
t tl i g +426 s, and PC data taken later on. The XRT data were pro- 
cessed following standard procedures dCapalbi et al. 20050), by 
using the version 0. 1 1 .6 of the dedicated pipeline xrtpipeline. 
Grade filtering was applied by selecting the 0-2 and 0-12 ranges 
for the WT mode and PC mode data, respectively. During the 
first orbit, the PC mode data show a count rate above 0.5 counts 
s" 1 thus likely a ffected by moderate pile-up. From point-spread 
function fitting dMoretti et al . 2005) we excluded a circular re- 
gion of 1 px radius from the extraction region of the source. 

The X-ray (0.3-10 keV) light curve from 87 s to 2.3 days 
after the trigger shows a steep initial decay up to 150 s after the 
trigger that is poorly approximated by a power law. Starting 
from 150 s after the trigger we find that the data are best fit- 
ted by a double broken power law model, with decay indexes 
a x ,i = 1.25 + 0.15 up to 1 ks, a x2 = 0.6 + 0.2 up to t b = 6 + 1 
ks and a x ,3 = 1.4 ±0.1 later on (x 2 /d.o. f. = 40/43). These val- 
ues are quite typic al for Swift/XRT afterglow light curves (e.g. 



Liang et al.l 120071) . We note that these results are marginally 
consisten t with the results f rom an automatic analysis of 318 
GRBs by lEvans et ail d2008h where, for GRB 0601 1 IB, a x ,i ~ 
1 from 175 s to 10 4 s after the trigger, corresponding to the 
average of our a X \ and a X 2- By fitting PC data only, a sim- 
ple power law and a broken power law model are rejected with 
> 99.99% and ~ 80% confidence levels, respectively. 

In order to investigate the spectral properties of each seg- 
ment individuated in the light curve, we extracted the 0.3-10 
keV spectrum of the X-ray afterglow at 4 different epochs, 
indicated with letters A, B, C, and D as plotted in Figure [3] 
Assuming an absorbed power-law model we found no evidence 
of spectral variation during the evolution of the afterglow, with 
a photon index of ~ 2 and an equivalent hydrogen column den- 
sity of Nh ~ 2 x 10 21 cirT 2 in addition to the Galactic one 



(N H ,Gai = 7.65 x 10 20 cm" 2 from |KalberlaetalJ 2005). These 
results are summarized in Table |4] Our tempor al and spectral 



results are in agreement with results found by Racusin et al 
(20083). 



2.3. Swift/UVOT 

The Swift follow-u p with UVOT started 84 s after the trigger 
dPerri et alJl2006albl: lHollandll2006l) . We analyzed the UVOT 
images starting with the reduction of raw data using the stan- 
dard SwiftUVOT software reduction tasks to create the sky 
images. We performed the analysis with the release 2.7.2 
(HEADAS version 6.3.1) and the 20070711 version of the 
CALDB. 

Our analysis involved the steps listed below. We performed 
the "standard" aperture photometry of the object in the sky im- 
ages for each single exposure and, as comparison, in the im- 
ages obtained as the sum of all the exposures in each filter. 
Photometry was performed using the UVOT-HEADAS tasks 
"uvotsource" and "uvotmaghist". For this purpose we first se- 



1 http://swift.gsfc.nasa.gov/docs/swift/analysis/ 
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Table 1. Swift/UVOT single exposure frames analyzed following standard procedure and photometry results. 



DATE 


Filter 


Exposure 


mag" 


mag 11 


K 


Log frequency 


Log (v Fy ) 


UT 




(s) 




(dereddened) 


(10- 28 ergcm-V'Hz-') 






2006-01-11 20:17:06 


V 


196.6 


17.59+0.11 


17.24 ±0.11 


46.19 ±4.67 


14.74 


-11.59 ±0.10 


2006-01-11 20:29:35 


V 


196.6 


> 18.86 


>18.51 


> 14.34 


14.74 


-12.10 ±0.00 


2006-01-11 21:52:36 


V 


118.8 


> 18.86 


>18.51 


> 14.34 


14.74 


-12.10 ±0.00 


2006-01-11 23:12:09 


V 


196.6 


> 19.29 


> 18.94 


>9.65 


14.74 


-12.27 ±0.00 


2006-01-11 20:20:33 


B 


196.6 


19.73 ±0.32 


19.27 ±0.32 


79.46 ± 2.34 


14.84 


-12.26 ±0.29 


2006-01-11 20:33:01 


B 


93.4 


> 19.24 


> 18.78 


> 12.48 


14.84 


-12.06 ±0.00 


2006-01-11 22:06:15 


B 


65.5 


> 18.98 


> 18.52 


> 15.85 


14.84 


-11.96 ±0.00 


2006-01-11 23:37:36 


B 


221.2 


> 19.84 


> 19.38 


> 71.80 


14.84 


-12.30 ±0.00 


2006-01-11 23:41:24 


B 


221.2 


> 19.67 


> 19.21 


> 83.97 


14.84 


-12.23 ±0.00 


2006-01-11 23:45:12 


B 


126.6 


> 19.43 


> 18.97 


> 10.47 


14.84 


-12.14 ±0.00 



a ] The upper limits are given at 3 cr confidence level. 



1999) with a value of E(B - V) = 0.11 mag dSchlegel et al 



19980 . The logarithm of the vF v calculated from this flux, i.e. 
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Fig. 3. The 0.3-10 keV light curve (blue data are in taken in 
Windowed Timing mode and red data are taken in Photon 
Counting mode by XRT). The continuum solid line shows the 
best fit temporal model (see §2.2). For each plotted temporal 
interval (A,B,C and D) we extracted an energy spectrum (Tab. 

a. 



lected the best source and background aperture regions from 
the first single exposure frame of the V filter. We applied the 
standard sizes for the source and background apertures, exclud- 
ing the nearby contaminating objects from the background ex- 
traction apertures. The results of the photometry of each single 
exposure in each filter is reported in Table 1 . We also report the 
dereddened flux obtained by correcting the observe d magnitude 
using a mean Galactic interstellar extinction law dFitzpatrick 



the Spectral Energy Distribution, is also reported. 

We further refined the analysis of the first single expo- 
sure frame of the UVOT image of filter V by applying "non- 
standard" apertures for the source extraction and by applying 
the aperture correction procedure in order to enhance the signal 
to noise ratio. Various "non standard" extraction regions, such 
as with larger radii or different shapes (for example an ellipti- 
cal region), have been tested, but generally give results similar 
to those obtained using the standard extraction radius. For the 
aperture correction procedure a first photometry on apertures 
of 1.9" (the FWHM of the V filter) was performed. Then 6 ob- 
jects not saturated in the field were selected and photometry 
was performed using FWHM apertures and standard apertures 
(5" radius) for each of them. The difference of the magnitudes 
in the two different apertures was evaluated for the stars se- 
lected and finally the weighted mean among all the differences 
was calculated. This value was subtracted to the GRB mag- 
nitude obtained with FWHM aperture. The final value of the 
GRB magnitude thus obtained integrating over all the exposure 
time is V = 17.55 + 0.08 mag that is consistent with our stan- 
dard analysis. The upper panel of Figure H] shows the field of 
the GRB in filter V, during the first single exposure frame of 
sky image. 

In order to increase the temporal resolution, the event data 
in all filters were also analyzed to extract the UVOT lightcurve 
of the object. Event data were first screened to filter out un- 
wanted contaminations using standard selection criteria, as ex- 
cluding the South-Atlantic Geomagnetic Anomaly and possi- 
ble anormal attitude variations, keeping the angle between the 
pointing and the Earth limb greater than 10 degrees and the 
same angle with "bright" Earth limb greater than 20 degrees. 
Moreover, only events classified as "good" were accepted. The 
lightcurve for each filter was extracted using the UVOT task 
"uvotevtlc", with different binning for the V filter (30 s bin) as 
compared to the remaining (70 s bin except for the White for 
which 100 s bin has been used). The photometry is calibrated 
using the photometric zero points in the Swift/UVOT calibra- 
tion data base. These zero point have a systematic error of 0.1 
mag that is not included in the values quoted in Table 2. 
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mid 


T exp (s) 


Filter 


mag b 


Amag b 


ref. 


29.65 s 


3.3 


R 


13.75 


0.06 


TAROT (1) 


32.95 s 


3.3 


R 


13.9 


0.09 


TAROT 


36.3 s 


3.4 


R 


14.05 


0.09 


TAROT 


39.65 s 


3.3 


R 


14.43 


0.13 


TAROT 


42.95 s 


3.3 


R 


14.63 


0.13 


TAROT 


46.3 s 


3.4 


R 


14.94 


0.12 


TAROT 


49.65 s 


3.3 


R 


15.11 


0.2 


TAROT 


52.95 s 


3.3 


R 


15.4 


0.14 


TAROT 


56.3 s 


3.4 


R 


15.58 


0.17 


TAROT 


59.65 s 


3.3 


R 


15.4 


0.22 


TAROT 


67.95 s 


13.3 


R 


15.79 


0.10 


TAROT 


81.25 s 


13.3 


R 


16.31 


0.13 


TAROT 


109.6 s 


30 


R 


16.46 


0.15 


TAROT 


146.5 s 


30 


R 


16.89 


0.20 


TAROT 


183.3 s 


30 


R 


17.53 


0.30 


TAROT 


85.5 s 


30 


V 


16.8 


0.4 


UVOT 


115.5 s 


30 


V 


17.2 


0.5 


UVOT 


145.5 s 


30 


V 


17.7 


0.7 


UVOT 


175.5 s 


30 


V 


17.6 


0.7 


UVOT 


388 s 


197 


B 


19.73 


0.32 


UVOT 


792 s 


300 


R 


18.9 


0.4 


MITSuME (2) 


792 s 


300 


I 


18.3 


0.4 


MITSuME <2) 


U.Zo days 


jyuu 


r> 
A 


Zj.Z 


U.J 


QTT 1 C/-1 

Kl 1-ljU 


12.19 days 


1440 


R 


>21.7 




Maidanak-1.5 


24.18 days 


2700 


R 


>22.3 




Maidanak-1.5 


27.20 days 


1800 


R 


>22.7 




Maidanak-1.5 


31.18 days 


3420 


R 


>22.5 




Maidanak-1.5 


16.29 days 


3360 


R 


25.1 


0.5 


ZTSh 


108.0 days 


1200 


V 


>21.6 




ZTSh 


26.3 days 


3600 


R 


>21.0 




OHP-T120 


37.1 days 


2700 


R 


>22.5 




NOT 


159 days 


1840 


R 


24.8 


0.4 


CFHT 


2.3 years 


1800 


H 


> 23.0 




CFHT 


5.3 months 


4800 




3800-8500A 


spectroscopy 


Keck/LRIS 



11 T mi d is the mean epoch of the exposure of duration T exp 

b Errors are given at l<x. The upper limits are given at 3 a confidence level. 

References. (1) all TAROT data are from lKlotz etail J2006h ; (2) lYanagisawa et al.l (2006) 

der to extract the maximum information on the OT, we care- 
fully analyzed the sum and median stacks of 13 images of 300s 
taken at RTT150 between 5.92 and 7.20 hours after the trigger 
(see also iKhamitov et al.l 2006). The afterglow is marginally 
detected with a signal to noise ratio of S/N=2. A PSF fit of the 
theoretical location of the afterglow gives R = 23.2 + 0.5. 

Further observations were performed with AZT-22 
(Maidanak Astronomical Observatory, Uzbekistan) 2.6m 
Shajn reflector (ZTSh, Crimean Astrophysical Observatory, 
Ukraina), OHP T120 (Observatoire de Haute-Provence, 
France) and NOT with only one detection at ZTSh 16.3 days 
after the trigger (see Table[2jl. 

Late time observations of the host galaxy are discussed in 
section 14. II 



2.4. Ground based optical follow-up 

Very early optical emission, measured with unprecedented 
temporal resolution of a few seconds by the TAROT ro botic 



observatory, has been presented by Klot z et al.l (12006). We 



quote here their results in Table [2] ROTSE-IIId observations 
that started 33 s after the trigger and detected a bright op- 
tical flare which was initially at ~ 13 mag and f ading as 



t ( " ±u u) , in agreement with TAROT observations (lYost et al 



20061 lRvkoffetalJl2009l) . The field of GRB 0601 1 IB was also 



observed with the three-color MITSuME 50 cm telescope at 
Okayama, Jap an starting at f fr ,> + 10.6 minutes in the Rc and Ic 
band images (lYanagisawa et al ] l2006h . 



Few hours after the trigger, GRB 06011 IB was ob- 
served with the Russian-Turkish 1.5-m telescope (RTT150, 
Bakyrlytepe, TUBITAK National Observatory, Turkey). In or- 
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Fig. 4. Top panel: R -band isophote curves of the host galaxy 
and the nearby bright galaxy from the CFHT image taken 159 
days after the trigger. Superposition of the UVOT image taken 
123 s after the trigger shows that the afterglow coincides with 
the host galaxy candidate from the CFHT image. Bottom panel: 
CFHT gray scale image is superposed to the UVOT isophotes 
of the afterglow. 

3. Multiwavelength analysis 

3.1. Temporal evolution 

The prompt emission of GRB 0601 1 IB was observed simulta- 
neously in the optical band ( R ) and in the gamma-rays (15-150 
keV), starting from 28 s up to 60 s (T90) after the trigger time. 
During this interval, the optical emission was found to decay 
steadily, in contrast with the gamma-ray emissio n which ex 



hibited flux variability of 2 orders of magnitude ( Kl otz et al 
2006h . 

The prompt optical emission from t tr ; g + 28 s to t tr i g + 200 s 
does not track the X-ray (0.3-10 keV) emission either. The X- 



ray light curve, from 



ttrig +87 s to 



Hng 



•150 s, shows a much 



steeper decay. However, starting from 200 s after the trigger the 
X-ray and optical decay indexes are consistent (see Tab.|4]and 
Fig. EJ. Indeed, the TAROT data, the MITS 11ME R -band ob- 
servations (taken 13 minutes after the trigger, Yanagisawa et al.1 
2006) and the RTT150 observations, show that the R -band light 
curve from 28 s to 7.2 hrs after the trigger can be modeled with 
a broken power law with best fit decay indexes qr 1 = 2.4 ± 0.2 
andffs, = 1.18+0.05 with abreakat ~ 75 s. The V -band decay 



Fig. 5. Temporal evolution of the prompt and afterglow flux at 
several wavelengths. The R -band early data are clearly not cor- 
related with the BAT prompt emission. The dashed curve is the 
R -band data best fit model while the solid curve is the X-rays 
best fit model. The dotted curve is the best fit X-ray model plot- 
ted over the R -band data. The three dotted vertical lines are the 
epochs at which a SED was extracted. 



index between ~ t lrig + 85 and ~ t lrig + 200 s is ay — 1.2 + 0.4, 
consistent with the /?-band decay index. The ZTSh detection 
(taken 16.3 days after the trigger) slightly exceeds the power- 
law extrapolation, likely due to the superposition of the under- 
lying host galaxy (see §4.1). 

GRB 0601 1 IB represents one of the few cases where color 
information during the early afterglow is available. The sim- 
ilarity of the V and R band decay indexes (within uncertain- 
ties) between 85 s and 200 s after the trigger is indicative of 
achromatic evolution of the late afterglow. However the large 
uncertainty affecting the V-band decay index ay and the prox- 
imity of R and V pass-bands prevent from more detailed inves- 
tigation. The lack of a continuous coverage in the optical band 
after about t tr i g +l ks, makes possible a consistency of the R - 
band light curve at late times with the X-ray shallow-to-normal 
behavior with a temporal break at about 6 ks after the trigger. 
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3.2. Spectral Energy Distribution 

In order to investigate on the origin of the prompt optical emis- 
sion, we analyzed the Spectral Energy Distribution (SED) from 
simultaneous BAT and TAROT observations. Two SEDs were 
extracted at 39 s after the trigger (from 28 s to 50 s, imme- 
diately before the second high energy prompt-emission peak) 
and at 55 s (from 50 s to 60 s, during the second peak). In 
the first interval, the available statistics from BAT data is too 
low to obtain spectral information and we could compute only 
the observed average 15-150 keV count rate, that is 1.6 x 10~~ 3 
counts s _1 . We converted this count rate into flux density as- 
suming the 790-averaged best fit spectral model (a simple 
power law with a photon index of T y = B y + 1 = 0.9 + 0.1 
see §2.1). Extrapolating the power law model back to low en- 
ergies, the TAROT flux density is about 3 orders of magnitude 
above the extrapolation of the gamma-ray flux density, with 
Fy/F opt ~ 10~ 3 and B 7t0pt ~ 0.6 (Fig. |6j. During the second 
peak of the burst (i.e. 55 s after the trigger), the available statis- 
tics enabled to extract a high energy spectrum best fitted by a 
simple power law with photon index of r r - B y + \ - 1.3+0.1. 
Fixing this value within the uncertainties, we found that the R 
-band extrapolation is consistent with the observed flux (Fig. 
[§). The FylF opt is now ~ 0.1 and B 7 ^ opt ~ B y = 0.3. 

These results assume no additional extinction beyond that 
corrected for in our Galaxy. However, from X-ray spectral anal- 
ysis there are hints of non-null absorption from the host galaxy. 
For example, assuming a Milky Way ISM, the rest frame A^ ost 
is exp ected to be jV ff / 1 . 8 X 1 2 1 mag cirT 22 jPredehl & Schmitt 
19951) . that for GRB 0601 1 IB is of about 1 mag or more since 
our Nh measure (Nh ~ 2 x 10 22 cm 22 assuming z = 0, 
§2.2) should be considered as a lower limit given the unknown 
redshift. We thus attempted to estimate the host galaxy dust 
extinction affecting the TAROT data by exploiting our mul- 
tiwavelength data set. We used Swift/XRT 0.3-10 keV and 
Swift/UVOT fi-band and V-band data as well as the TAROT 
/?-band and MITSuME /-band data. Since we have only one 
detection in the fi-filter, the most sensitive wavelength to dust 
extinction, we decided to extract the SED at the time of the 
fi-filter detection, that is at 388 s after the trigger. /-, R- and 
V-band data have been extrapolated at f, ns; +388 s assuming the 
/?-band decay slope at that time, that is 1.18 + 0.05, and then 
corrected for Galactic extinction, that in each filter is A#=0.48, 
Ay=0.37, Ar=0.30 and A/=0.22 mag, using a standard galactic 
extinction curve. 

The optical spectral index is unusually steep (B opl > 2). For 
comparison, other optical spectral indexes of GRBs with early 
optical steep-to-shallow light curve for which spectral data are 
available during the shallow pha se are e.g.B op , = 0.17 for GRB 
0803 19B JWozniak et alJ Eooi and B opt = 0.95 GRB 061126 



dPerley et al.ll2008|) . This lends additional support to our infer- 
ence of significant extinction in the host galaxy. 

The consistency of the decay slopes in optical and X-rays 
at f,rig+388 s suggests a common cooling regime of the elec- 
tron populations re sponsible for the observed emission in the 
two energy ranges (ISari et al.ll 1998b . Specifically, the flux de- 
cay rate and the X-ray spectral index are consistent with the 
fireball expected closure relationship for an electron popula- 



tion in slow cooling regime and emitting at frequencies above 
the cooling frequency v > v c . The expected spectral and de- 
cay indexes are p/2 and (2 - 3p)/4 respectively, where p is the 
electron distribution spectral index. We find a good agreement 
with the measured indexes for p ~ 2.3. The optical and X-ray 
SEDs, both extracted at f„ Ji; +388 s, were therefore fitted simul- 
taneously assuming a simple power law model with spectral 
index fixed (within 1 cr) at the value obtained from X-ray data 
analysis only, that is ax = 1.2 + 0.1 (e.g. interval B in Table|4]l. 
To account for dust extinction from the host galaxy, we tested 
three different extinction curves. 

Assuming d ifferen t redshifts, the best fi t among Galactic 
lCardelli et all Il989h . SMC H Il992h and Starburst 



Calzetti et al.l 1 19941) extinction curves was obtained with 
the latter at z = 2 with rest frame visual extinction A r v ~ 2 mag 
(Fig. El Tab. |5J. We have assumed a ratio of total to selective 
extinction of Ry = Ay/E(B - V) = 3.1 for the Galactic and 
the Starburst extinctio n curv es while Ry - 2.93 for the SMC 
extinction curve (Pel 1992). At z < 2 acceptable fits were 
obtained only with the Starburst extinction curve while at z > 2 
comparable good fits were obtained also assuming the Galactic 
extinction curve up to z - 2.5 - 2.7, above which the Lyman 
a cut-off due to inter-galactic neutral hydrogen absorption 
would provide a B-band flux density more than 2<x below the 
observed one. 

At z — 2, the observed /?-band corresponds to the UV band 
in the rest frame (central wavelength ~ 0.25//m), with best fit 
rest frame extinction of Ajjy - 2.5 + 0.7 mag. We note that at 
z — 1 and z — 3 the corresponding rest frame UV extinction is 
2.7 and 1.7 mag, respectively. Indeed, if at z > 2, part of the B- 
band attenuation may be due to Lyman a cut-off and not only 
to dust extinction. Despite the redshift uncertainties, these find- 
ings indicate that the observed optical flux is attenuated by host 
dust extinction. Correcting the observed R -band magnitude of 
~ 2.5 mag, that is the best fit value with a starburst extinction 
curve, the optical flux density increases by about a factor of 10. 
The above analysis of the gamma-to-optical SED would now 
provide F 7 /F opt ~ 10~ 4 and B %opt = 0.80 + 0.05 at 39 s after 
the trigger and F 7 /F opt ~ 0.01 and B %op , = 0.50 + 0.05 at 55 s 
after the trigger. 

These results make GRB 0601 1 IB consistent with the few 
GRBs with prompt optical detection for which the extrapola- 
tion of the prompt gamma-ray flux to the optical bands assum- 
ing a simple po wer law un der-predicts the observed flux (cases 
'B' and 'D' in I Yost etaD 20071. Together with their different 
evolution in time, these findings suggest two separate origins 
for the prompt gamma-ray and the optical emission compo- 
nents.^ 

In lKlotzet al.N2006l) a value of A r v ~ 4 mag at z ~ 1 was 



reported from our XRT and TAROT data at ~ t tr i g + 140 (as- 
suming a simple underlying power-law spectral continuum and 
a Galactic extinction curve). This is significantly higher than 
our estimate at f /r , ? +388s, possibly indicating an active dust 
destruction mechanism or the presence of a spectral break be- 
tween the two bands at that time. In the latter case, the decaying 
flux behavior and the similarity of the decay rates observed at 
optical wavelengths and in X-rays from ~ 200 s after the trig- 
ger suggest that the synchrotron cooling frequency crossed the 
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rest frame UV band between t trig + 140 s and t trig + 200 s after 
the trigger. 
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Fig. 6. Spectral Energy Distribution at ttrig + 39 s and t tr i g + 55 
s after the trigger wh ere optical data hav e been corrected 



for Galactic extinction (Schlegel et al.lll99 8) but not for host 



galaxy extinction and should be considered as lower limits. The 
solid curve represents the optical-to-gamma-ray slope while 
the dashed lines individuate the best fit power law model from 
gamma-ray data only within 1 cr (at t tr i g + 39 s we have assumed 
the Tgo integrated spectrum best fit model). 



4. The host galaxy of GRB 0601 1 1 B 

4.1. Imaging 

Two deep exposures have been taken with the Canada-France- 
Hawaii Telescope (CFHT0) 5.3 months and 2.3 years after the 
burst with the R and H band filters, respectively. We detect a 



2 Based on observations from director discretionary time obtained 
at the Canada-France-Hawaii Telescope (CFHT) which is operated 
by the National Research Council of Canada, the Institut National 
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Fig. 7. Spectral Energy Distribution at f /r „, + 388 s where optical 



data have been corrected for Galactic extinction ( Schlegel et al 
119981) . The solid curve is the optic al-to-X-ray best fit model ob- 
tained at z = 2 assuming a Starburst extinction curve (see §3.2) 
and the open red circles are the expected flux density values in 
each energy band. 

Table 3. Best fit parameters obtained for the Spectral Energy 
Distribution extracted at t tr i gg + 388s assuming a simple power 
law model. 



Redshift 


Extinction Model 


mag 


x 2 /y 





Galactic 


5.2 + 0.2 


19.8/7 




SMC 


5.5 + 0.2 


20.7/7 




Starburst 


4.1 ±0.3 


5.2/7 


1 


Galactic 


2.2 ±0.2 


14.0/7 




SMC 


2.1 ±0.2 


15.2/7 




Starburst 


2.8 ±0.2 


3.5/7 


2 


Galactic 


1.5 ±0.1 


5.1/7 




SMC 


2.1 ±0.1 


15.2/7 




Starburst 


2.3 ± 0.2 


2.8/7 



persistent source in the /?-band, consistent with the afterglow 
position of GRB 06011 IB, with magnitude R = 24.8 + 0.4. 
The host is 6.3" off-axis from a nearby, bright spiral galaxy. 
From the positions derived from the relative astrometry per- 
formed on both the UVOT and CFHT images, we found that 
the afterglow and the host galaxy have an offset of less than 
0.2 arcseconds (Fig.©. The GRB 06011 IB field was also ob- 
served in the //-band for an exposure time of 30 minutes with 
the WIRCAM wide-field NIR camera at the CFHT on May 18, 
2008 (2.3 years after the burst). Data were preprocessed and 
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Epochs 


Instrument 


Energy range 


Decay Index 


Spectral Index 


T - T, rig (s) 










28-50 


BAT 


15-150 keV 






50-60 


BAT 


15-150 keV 




0.3 ± 0.2 


28-75 


TAROT 


R -band 


2.4 ± 0.2 




75-26000 


TAROT+MITSUME+RTT 


R -band 


1.18 ±0.05 




85-200 


UVOT 


V -band 


1.2 ±0.4 




A(87-150 s) 


XRT/WT 


0.3-10 keV 


3.8 ±0.5 


1.0 ±0.2 


B(150- 1000s) 


XRT/WT+PC 


0.3-10 keV 


1.25 ±0.15 


1.2 ±0.2 


C(4350-6800s) 


XRT/PC 


0.3-10 keV 


0.6 ± 0.2 


1.3 ±0.2 


D(9930- 185650s) 


XRT/PC 


0.3-10 keV 


1.4 ±0.1 


1.1 ±0.1 



calibrated by "I'iwi" pipeline at CFHT, and stacked with the 
SWARP software. We could not detect any source at the loca- 
tion of the burst down to the limiting magnitude of 23.0 (22.0) 
for a signal to noise ratio of 1.5 (3.0) (AB magnitude system). 

From the CFHT image and considering that at the epoch 
of its exposure the afterglow has definitively faded away, it is 
clear that GRB 0601 1 IB occurred within a faint galaxy. 

4.2. Spectroscopy 

A spectroscopic observation of the host galaxy has been per- 
formed the by LRIS instrument on the Keck I telescope on 2006 
July 26 (UT). The 600/8500 grism and 600/4000 grating with 
the 560 dichroic were used in order to allow the covering of a 
wide wavelength domain (between 3800 and 8500 A). The slit 
width was fixed to 1.5 arcsec and it was oriented in order to 
cover both the GRB host galaxy and the bright nearby galaxy 
(see §4.3). The total integration time was 2 x 2400 s. Due to 
the position of the object in the sky, the airmass was relatively 
large during the observation, about 1.6. Data were reduced us- 
ing the LRIS rutines within the IRAF package. The two arms, 
Blue and Red, were reduced separately. In particular, the over- 
scan region was removed and the bias was subtracted by using 
the Irisbias task, which manages the dual amplifier (two-amp) 
readout mode. Since the two LRIS arms do not have exactly 
the same spectral resolution, the two spectra were therefore re- 
binned to the same wavelength resolution of 1 A and finally 
combined together to produce the final spectrum. The junction 
region between the two arms, around 5800 A, presents anoma- 
lous features that were cleaned 'manually' (Fig. [8]). 

No clear emission or absorption lines are visible in the 
spectrum. No feature is compatible with a galaxy located close 
to, or at the same redshift than, the nearby large galaxy (see be- 
low). By applying a simple boxcar filter over 15 pixels, a weak 
continuum is observed, which is compatible with the one from 



Rest frame wavelength (A) 
1200 1400 1600 1800 2000 2200 2400 2600 2800 3000 



a sta rburst spiral with 0.39 < E(B - V) < 0.50 dKinnev et al 
[1996J) at a redshift z ~ 1.9. 

4.3. The bright nearby galaxy 

A much brighter galaxy is present in our imaging as well, cen- 
tered at an angular offset of only 4.5 arcseconds from the GRB 
host galaxy. Astrometry based on CrAO R-band observations 
taken 16 days after burst and USNO-A2.0 reduction located 
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Fig. 8. Upper panel: optical spectrum of the host. The red curve 
is the template model. Bottom panel: optical spectrum of the 
nearby galaxy taken with the Keck Telescope. 



the bright galaxy at coordinates (J2000) R.A.=19'' 05 m 43.06*, 
Dec.=+70d 22' 29.7" with uncertainty of 0.2" (Fig.HJi. We per- 
formed first photometric and then spectroscopic observations 
of this galaxy in order to check whether any connection with 
the host (e.g. a satellite galaxy) and/or the burst itself (e.g. a 
halo burst as if it originated from an old binary system) could 
exist. Photometric observations have been performed at the 
NOT telescope about four and eight months after the burst. U, 
B, V, R, I and J magnitude are quoted in Table [5] T he best pho 



tometric redshift obtained with HyperZ code ( Bol zonella et al 
20001) is z — 0.22^^. Its type is compatible with an early spi- 
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Table 5. Photometry of the nearby galaxy (not corrected for 
Galactic extinction). 



Filter 


magnitude 


Amag 


Exposure 

(s) 


U 


23.34 


0.3 


1500 


B 


22.54 


0.06 


810 


V 


21.02 


0.03 


360 


R 


20.24 


0.02 


360 


I a 


19.70 


0.03 


180 


J 


17.63 


0.05 


2880 



Filter I„„ is centered at J=7970A(FWHM= 
A) #12 interference I filter more details 
http://www.not.iac.es/instruments/filters/filterlist.html 



157 
in 



ral galaxy, Sb. From J band images taken at NOT with WFcam, 
we find marginal evidence of two spiral arms. The angular scale 
at the measured redshift is 3.76 kpc/arcsec, making the size of 
this galaxy close to 5.5" x 3.76=20.7 kpc. Hereafter, we as- 
sumed a flat cosmology with Hq = 71 km s Mpc, O a = 0.73, 
Cl m = 0.27 (Wright 2006). From spectroscopic analysis per- 
formed of the Keck spectrum (see previous section), the galaxy 
displays several emission and absorption lines (Fig. [8}, allow- 
ing redshift measurement of z = 0.239 (1 % accuracy). 

4.4. Limits on the SN component 

No detection of any associated supernova (SN) down to R = 
25.1 was observed in the ZTSh telescope observations taken 
16 days after the GRB (Tab.©. Assuming a s GRB-SNe tem- 
plates the cases of SN 19941 (GRB 2121 1. iDella Valle etal 



2004) and SN1998bw (GRB 980425. iGalama et al.l 1998) and 
considering a host galaxy dust extinction in the V band of about 
~ 1 the lack of any SN detection in our images implies a red- 
shift lower limit of about 0.4 for a SN1994I-like event and of 
about 0.7 for a SN1998bw-like one. 

4.5. Redshift estimate 

In this section we attempt to estimate the distance of GRB 
06011 IB following some considerations based on our multi- 
band observations of both the host and the afterglow. 

The host galaxy spectrum shows a weak continuum that is 
comp atible with a starb urst spiral with 0.39 < E(B - V) < 
0.50 dKinnev et al.lll996h at a redshift ~ 1.9. Unfortunately, no 
emission or absorption lines are present (Fig. [8j. The two most 
common lines observed in long GRB host galaxies are the [Oil] 
line and the Ha. Although the lack of these lines cannot give 
an absolute restriction to the redshift range, the lack of [Oil] (A 
3727 A) emission line may imply z > 1.5 while the lack of 
Ha(A 6563 A) emission line may imply z > 0.4. 

The best-fit redshift from the afterglow SED was z — 2, 
consistent with the result obtained from fitting the host spectral 
continuum. We also used the X- ray light curve to est imate the 
distance to this burst: following iGendre et al] (I2008al) we esti- 



mate z = 2.2 + 1.0, consistent with the above results. Finally, 
comparing the host galaxy apparent magnitude of R ~ 25 with 
the typical magnitudes observed for a sample of long GRBs 
hosts, we find that it is likely at z ~ 1-2 (cf. for instance 
Figure 9 of Berger at al. 2007), again consistent with previous 
findings. 

We further investigate on the distance scale of this burst 
by comparing the energetics and the host galaxy properties ob- 
tained assuming different redshifts, with the typical values ob- 
served for long GRBs. In Table|6]we have computed the size of 
the host (assuming an angular size of 1", corresponding to the 
seeing of CFHT observation) and its absolute magnitude and 
the burst energetics at several different redshifts. The burst en- 
ergetics was estimated with the isotropic equivalent energy Ei SO 
that we have evaluated through the observed 100-1000 keV flu- 
enceS = (5.6+0.8)10~ 6 erg cirT 2 (see §2.1) as S x4nD 2 L l{\+z) 
where Di is the luminosity distance. The absolute magnitudes 
have been computed from the observed R magnitude, consider- 
ing that for z > 0.5 it samples the rest frame B filter magnitude 
while for z > 1 we assumed U - B ~ 0.75 mag as typical U — B 
color of starburst galaxies (e.g. lCoil et al . 20071 |Berger||2009). 



Table 6. Host galaxy and burst properties at different redshifts. 
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d a 


M\ 


F 


E peak J 






kpc 


mag 


10 52 erg 


MeV 


deg 


0.24 


3.76 


-15.4 


0.07 + 0.01 


0.6 ± 0.3 


4.5 


1 


8.0 


-18.4 


1.5+0.2 


0.9 ± 0.5 


2.6 


2 


8.5 


-20.3 


5.8 + 0.8 


1.4 ±0.8 


1.9 


3 


7.8 


-21.5 


12 + 2 


1.8 ± 1.0 


1.5 



a host size 

b absolute magnitude of the host 

c jet opening angle assuming that the temporal break observed in 
the X-ray light curve at a~ 6 ks has a jet origin 



Long GRB host luminosities cover a wide range of val- 
ues, most of them between Mb = - 17 and Mb = -22 mag, 
with a peak at about Mb ~ -18 mag (Berg er et alj|2007l) whil e 
the typical host sizes is ~ 1.7 kpc dWainwright et al] l2007h . 
Considering that the host angular sizes quoted in Table|6]should 
be considered as an upper limit due to the seeing at CFHT at 
the epoch of the observations, the host size of GRB 06011 IB 
is consistent with typical long GRBs host sizes at all the tested 
redshifts while its absolute magnitudes would be too faint at 
the redshift of the nearby galaxy. From burst energetics con- 
siderations, Table shows that at z = 0.24, GRB 06011 IB 
would be an outlier in the Amati E pea ^i versus Ei so correlation 
(Fig. |5). We note however that even at higher redshift GRB 
06011 IB seems to deviate from the correlatiorH Other long 
GRBs de tected by BATSE are inconsistent with the A mati cor- 
relation dNakar & Piranll2005l: iBand & Preecell2005l) and GRB 
06011 IB may fall into this category (Fig. |9). However, the 
large uncertainties affecting the peak energy measure prevented 
us to reach any firm conclusion. 



3 This may be the reason why the pseudo-redshift (Atteia 2004) 
computations provides the uncomfortably high value of pz ~ 8.4. 
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Fig. 9. GRB 0601 1 IB is plotted (in red) in the E peak vsE iso plot 
at four different redshifts: z = 0.24,1,2 and 3, the latter being 
the maximum redshift compatible with our observations (see 
§5.2). Other GRBs data and the best fi t correlation whithin 2 cr 



are taken are from lAmati e t al. (2008) 



All of these considerations argue against any relation with 
the nearby, bright galaxy. Consistent with this conclusion is the 
lack of any supernova detection that imposes a redshift larger 
than z ~ 0.4 - 0.7. For the latter, we assumed an observed mag- 
nitude of R>25.1 mag for the SN; that is, the same brightness 
as the source detected in the ZTSh telescope images taken 16 
days after the burst (see §4.4). 

An upper limit on the redshift can be obtained from the lack 
of evidence of Lyman alpha cut-off due to the absorption from 
neutral hydrogen in the intergalactic medium in the spectral 
continuum of the host galaxy and by the B-band detection, 
providing z < 3. 

Summarizing the above considerations, there is evidence 
against a possible relation of GRB 06011 IB with the nearby, 
bright galaxy and independent estimates converge to z ~ 1-2. 
Having GRB 0601 1 IB at this redshift would result in having a 
standard long GRB host galaxy. 

4.6. Outflow collimation 

Using our redshift estimates above, we can compute the pos- 
sible collimation of the outflow Gj assuming a jet origin of the 
temporal break tj, observed in the XRT light curve at about 6 
ks. The angle estimations computed assuming a constant ISM 
density of n ~ 5 cirT 3 and radiative efficiency ~ 0.2 are pre- 



Fig. 10. Early optical emission light curves as if all bursts 
were at z = 1 of GRB 06011 IB (filled circles), GRB 990123 
(open s quares. lAkerlofetal.1 199 91. GRB 021211 (open tri- 



angles, iFox et al] 20 03. iLi et alJ 2003), GRB 060117 (open 
circles, Ijelfnek et al] 2006), GRB 061126 (open pentagons, 



Perlev et al. 2008), GRB 0803 19B (dashed line, best fit from 



Racusinetal. 2008b) and GRB 090102 (cyan pentagons, 
Klotz et all 2009). For these GRBs the Reverse Shock model 
has been claimed to explain the prompt optical afterglow (see 
Tabled}. 

sented in Table 5. It occurs that the jet opening angle is within 
the distribution of angles o f sample GRBs with well defined jet 
break in X-ray light curve dRacusin et al . 2008a). Alternatively, 
if the jet break epoch is beyond the end of XRT observations, 
that is tj, > 200 ks, then 6j > 6 deg (z < 2) which is marginally 
consistent with the distribution. Thus we can support the break 
at tb ~ 6 ks as a real jet break. This result is in line with 
the anticorrelation found betwee n the X-ray aftergl ow shallow 
phase durations and the redshift dStratta et al . 2009), that, if in- 
terpreted as a selection effect, may reflects an anticorrelation of 
the X-ray shallow phase duration with the burst energy. Indeed, 
the latter is expected if the end of the shallow phase originates 
from a collimated outflow. The burst energy corrected for col- 
limation at z = 1-2 would be E y ~ 1 - 2 x 10 49 erg. 

5. The nature of the prompt optical emission 

The nature of the prompt optical emission of GR B 06011 IB 



has b e en a l ready investigated in past works (e.g. iKlotz et al 



200o| |WeJ 120071) . A reverse shock (RS, e.g. Meszaros & 
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Rees 1994, 1 997) e mission has been invoked for this burst by 
Klotz et al.1 (120061) . primarily based on the TAROT observa- 
tions of very steep early-time optical decay (a = -(2.4 + 0.1) 
in R-band, see §3.1). Further evidence includes the observed 
lack of temporal correlation between the simultaneous BAT and 
TAROT emissions and the similarities of its optical light curve 
with GRB 990123 and GRB 02121 1 (for which a reverse shock 
origin has been proposed to explain their optical light curves). 
A decay index of ~ -2 has generally been considered as the 
evidence of a RS origin. Some concerns however have been 
addressed to the unknown RS zero-time that, if different than 
the trigger time, would make the real decay rate of the emission 
different from tho se measured, leading to different conclusions 
( Panaitesc"ull2009l) . A simultaneous spectral and temporal anal- 
ysis would be the best way to investigate the nature of the ob- 
served optical radiation. Unfortunately, for GRB 06011 IB we 
do not have any spectral information in the optical regime dur- 
ing the steep decay. In this work we further investigate on the 
nature of the prompt optical emission by exploring its optical- 
to-gamma ray spectral energy distribution. We also compare 
GRB 06011 IB with other GRBs showing early optical emis- 
sion with similar temporal properties. We then attempt to con- 
strain the physical parameters of this burst using our combined 
early-through-late-time observations, including our constraints 
on the redshift and on the host-galaxy dust extinction. 

Thanks to our multi-band analysis we were able to con- 
firm the presence of non-negligible extinction from dust in the 
burst environment. By correcting the TAROT data for the esti- 
mated amount of host galaxy dust extinction, we could evalu- 
ate the intrinsic optical-to-gamma-ray flux density ratio during 
the prompt emission. The obtained values are orders of mag- 
nitude above unity, quite unlike most GRBs with simultane- 
ous optica l counterpart which have been reported to date (see 
Yost et al.1 2007). Most importantly, the optical flux density is 
well above the extrapolation obtained from the best fit spec- 
tral model in the gamma-ray regime. This result, as well as the 
smooth, fast decay of the optical light curve simultaneous to 
the strong gamma-ray flux variability, enabled us to firmly ex- 
clude a common emission mechanism at high and low ener- 
gies. To compare GRB 0601 1 IB with other GRBs with "steep- 
to-shallow" decay observed during early optical emission, we 
plotted the optical magnitudes as a function of time as if all 
bursts were at z = 1 (Figure ITOb . To compute the expected 
magnitudes at z = 1 we simply apply the distance modulus 
definition without taking into account any k-correction and the 
host galaxy extinction correction since for almost all the plot- 
ted bursts we have no color information during the early op- 
tical emission and very uncertain or no host galaxy extinction 
measurements. For GRB 0803 19B, given the huge amount of 
optical data available, we decide d to plot only the we ll con- 
strained best fit light curve model dRacusin et al.ll2008al) in the 
period of interest for our study (that is, when the flux features 
the steep-to-shallow decay). Inspection of the figure demon- 
strates the strong similarities in the initial steep decay indexes 
(all within ~ 2 and ~ 2.5, see Tab. [7]). This evidence possi- 
bly indicates the presence of a common starting time of the 
observed prompt optical emissions that, if coincident with the 



trigger time, provides good consistency with the RS predicted 
decay slopes. 

In the context of the internal-external shock model (e.g. 
Sari & Piranlll997[) . the RS emission mechanism is a conse- 
quence of the external shock that the outflow produces with 
the surrounding interstellar matter, along with a forward shock 
(FS) that generates the afterglow emission. As for the FS, the 
RS accelerates electrons in the ou tflow into a power law energy 



distri bution with index p (e.g. ISari et al.lll998t ISari & Piran 
1 1999b . The light curve at a given frequency depends on the tem- 
poral evolution of the injection frequency v,„, the c ooling fre- 



quenc y v c and the synchrotron peak flux F VMax (e.g. lSari et al 



Il998l) . Kobayashi (2000) derived a simple analytic solution of 
the RS light curve in two limits: the thick shell limit, when the 
RS becomes relativistic crossing the outflow and the thin shell 
limit, when the RS remains Newtonian. From an observational 
point of view the thick shell limit predicts that the RS peaks 
during the high-energy burst while a RS peaking after the end 
of the burst is predicted from the thin shell limit case. The com- 
bina tions of the RS a nd FS light curves have been classified 
bv lZhang et al. I d2003h as of type 1,11 and III depending on RS 



and FS relative peak times and peak flux intensities. iKlotz et al 



(2006) pointed out that the optical prompt emission o f GRB 
0601 1 IB can be classified as type II (IZhang et al . 2003) where 
the RS emission is dominating over the peak of the FS compo- 
nent after the FS peak time. 

In order to compare observations with theoretical predic- 
tions, we first attempt to constrain the fireball initial Lorentz 
factor r from the t i me at which we o bserve the prompt optical 
decay, (ISarilll997t IZhang et al.ll2003l) . Since the optical light 
curve is already fading at t — 28 s after the trigger, the rest 
frame deceleration time is at f</ < f/(l + z); that is, the fireball 
has given the ISM an amount of energy comparable to its ini- 
tial energy at = [3E /32jmm p c 5 r & ] 1 ^ < f/(l + z). Assuming 
an interstellar matter density of n ~ 5cm 3 and a redshift of 
Z = 1 — 2 and estimating the energy of the shell E with Ej SO , we 
find F >260-360. This result is consistent with past estimates 
of the fireball Lorentz factor of about a few hundreds (e.g . 
Molinari et ap2007l: iGomboc et al.ll2008l: iFerrero et al.ll2008h . 

2006) GRB 06011 IB 



As already mentioned by Klotz et al 



falls in the thick shell limit since the optical prompt emission is 
already fading before the end of the prompt, high-energy burst. 
In this limit, the deceleration time should coincide with the 
time at which the RS crosses the shell (e.g. Ko bayashi 2000; 
Zhang et al 2003) that is estimated with the high-energy burst 
duration. We use 7go ~ 25 s observed in the Suzaku/WAM 
energy range (100 - 1000 keV) to estimate the burst duration. 
The injection frequency, the cooling frequency and the peak 
power at the shock crossing time fo r the thick shell limit are 
dSari & Piranlll999HKobavashill2000h : 



r ~73x10'V f 1/2 K 1/2 r 2 H? 



v c ~ 4.5 x 10 15 e B 3 ^ I £ 52 I/2 »5 1 (r 1 oo/(l +z))" 1/2 //z 



1/2, 



(1) 

(2) 



F VMax ~ \30D-^e]ll M E^n'l A T-^{T m) l{\ + zW^mJy (3) 

where Dl.2% - £>l/10 28 cm and £52 = £,, o /10 erg are the 
luminosity distance and the outflow energy computed at z = 2 
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(Tabled, e Bj o.oi = eg/0.01 and e e>06 = e f /0.6 are the fractions 
of the burst total energy that goes into the electrons and mag- 
netic field respectively, ng = n/5 cirT 3 is the interstellar matter 
density, Tioo = 7W100 s > an d F3oo = F/300 is the outflow 
Lorentz factor. Assuming the specified values, from equations 
(1) and (2) we find that at the RS crossing time the observed 
R-band frequency (v ~ 1.2 x 10 15 Hz in the rest frame), is be- 
tween v m (7.3 x 10 14 Hz) and v c (4.5 x 10 15 Hz). Equation (3) 
shows that the RS peak flux density is expected at ~130 mJy 
for z — 2. From Figure 5 we see that at the beginning of the 
prompt optical decay, the flux density is ~ 10 mJy. This value 
however is not corrected for the dust extinction in the host. The 
latter is highly uncertain but we find evidence that is not neg- 
ligible (see §3.2). From our optical-to-X-ray best-fit obtained 
assuming a starburst extinction curve, we measure a rest frame 
UV extinction of about 2.5 mag that makes the corrected mea- 
sured flux density consistent with the expected peak flux. This 
result indicates that the RS peak time should not be too far from 
the TAROT start observing time and that the outflow Lorentz 
factor is close to the estimated lower limit value. 

From (1) and (2) we obtain that v„, < v c at the RS crossing 
time. Therefore, for the TAROT observations we can assume 
a slow cooling regime. For v m < v < v c as observed, F v = 
F v Jv/v m y (p - l)/2 oc r 0. Since v m « r 73/4S ~ r 3/2 and F Vm oc 
1 dKobavashf 2000). we expect a = (3p + l)/4 that 



-47/48 



is consistent with our measured decay index if p — 2.9 + 0.3. 
We note that this value is somewhat higher than what derived 
from the FS component (see §3.2). Different values of p for the 
FS and the RS hav e been invoked as a possible scenario also 
for GRB 061 126 bv lPerlev et al] d2008l) . 

It has been argued that the RS emission mechanism can- 
not be the dominating mechanism to explain the bright optical 
early emission due to the lack of this c omponent in many other 
GRBs despite the expected high flux dRoming et alj |2005). We 
note that the prompt optical emission in GRB 06011 IB ends 
very early compared to the other events (Fig.fTOt. opening the 
possibility that this component may remain undetected in a sig- 
nificant fraction of GRBs. In addition, the magnetization of 
GRBs outflows pla ys a cruc ial role in suppressing the RS opti- 
cal emission (e.g. lFanll2008h . 

In a scenario alternative to the RS, the prompt optical emis- 
sion may be produced by internal shocks different from those 
responsible for the gamma-ray emission, which can be formed 
by collisions of different shells: these shells may be easily dis- 
rupted by other shells, thus ex plaining wh y bright optical emis- 
sion is not common in GRBs dWeill2007l) . For GRB 0601 1 IB, 
Weil (120071) computed that the rest frame cooling frequency 
should be at energies v c ~ 0.3 - 0.7 keV, that is between the 
X-ray and the optical energy domains, at z = 1 - 2, 28 s after 
the trigger and assuming a Lorentz factor of the merged shell of 
800. The v c expected value at the shell crossing time is consis- 
tent with our observations assuming that v c remains at energies 
higher than the optical up to f,„ g + 200.? after which the de- 
cay indexes in the two energy domains are consistent among 
themselves (see §3.1). However, the v c theoretical estimation 
depends on the highly uncertain fireball Lorentz factor to the 
power of eight. If for example we assume F = 300, the cooling 
frequency would had been well below the /?-band at t tr i g + 28.?. 



For this reason and for the numerous pieces of evidence sup- 
porting the RS emission exposed above, we do not favor the 
internal shock interpretation. 

Finally we mention that the RS emission mechanism has 
been proposed to solve the difficulties of the 'standard' for- 
ward shock model to explain the canonical (steep-flat-normal) 
evolution of the X -ray emission dUhm & Beloborodovll2007t 
Genet et alJl2007ah . GRB 06011 IB represents one of the few 
cases where both an opti cal flash and a 'canonical' X-ray light 
curve dNousek et al.l20 06) are observed, enabling to test the RS 
in both the energy ranges. For example, to reproduce the ob- 
served X-ray light curves a long-lived RS emission is required, 
which is difficult to reconcile with bright optical early emission 
dGenet et alj|2007al) . 



Table 7. Best-fit decay indexes of the early optical emission for 
all GRBs with well monitored optical counterpart that shows a 
steep-to-shallow behavior (Fig.fTOb. 



GRB 


Filter 


»i 


a 2 


' break 

(s) 


z 


990123 (1) 


V 


~ 2.5 


~ 1.5 


700 


1.60 


02121 1 (2) 


R 


1.97 ±0.05 


1.01 ±0.02 


300 


1.0 


060117 (3) 


R 


2.49 ± 0.05 


1.47 ±0.03 


300 


l a 


061126 (4) 


R 


2.0 ±0.3 


0.76 ± 0.06 


780 


1.1588 


06011 IB 


R 


2.4 ±0.2 


1.18 ±0.05 


75 


2 


080319B <5) 


R 


2.49 ± 0.09 


1.25 ±0.02 


900 


0.937 


090102 (6) 


R 


1.83 ±0.13 


1.1 ±0.3 


400 


1.547 



For GRB 0601 17 since the redshift is not known, we imposed z=l. 



References. (1) lAkerlof et alJdl999h: (2)lHolland et alj (120041): (3 ) 
IJelfnek etafl I2OO6I): (4 ) |Perlev et ail (120081): ( 5) lRacusin et ail 
( l2008iJ) : (6) lKlotzet al]te009b . 



6. Summary and conclusions 

The presence of a bright optical emission during the prompt 
phase makes GRB 0601 1 IB an extremely interesting burst and 
for this reason an extended multi-wavelength campaign has 
been performed in order to provide further constraints to the 
physical origin of the observed optical emission. Hereafter we 
summarize our main results and conclusions. 

Our late optical observations enabled us to identify the host 
galaxy of GRB 06011 IB. The featureless host spectrum, the 
non detection of any associated supernova 16 days after the 
burst as well as the inferred host luminosities at different red- 
shift, allowed us to put some constraints on the distance scale 
of GRB 0601 1 IB. Independent estimates converge to z ~ 1-2. 

From the afterglow SED we find that a non-negligible host 
galaxy dust extinction affects the observed flux in the optical 
regime. After correcting for this effect, the optical-to-gamma- 
ray flux ratio during the prompt strongly suggests a separate 
origin of the optical component fr om that one at high energies. 
This was already pointed out by dKlotz et alj|2006l) from the 
observed high energy prompt light curve not tracked by the si- 
multaneous optical one. The similarities of the measured opti- 
cal decay index of GRB 0601 1 IB with other observed prompt 
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optical emissions suggest a common zero-time origin of the 
emission mechanism that, if coincident with the burst trigger 
time, makes the measured decay indexes consistent with re- 
verse shock scenario. 

From the time at which we observe the TAROT light curve 
starting to decay, we estimated the fireball initial Lorentz factor 
dSarill 19971) as F >260-360 at redshift z ~ 1 - 2, in agreement 
with past results obtained for other bursts. Assuming F ~ 300, 
the observations are in good agreement with theoretical predic- 
tions from the reverse shock emission in the thick shell limit. In 
particular, we find evidence that the prompt optical emission is 
in the slow cooling regime with v m < v < v c and that the peak 
flux is very near or coincident with the observed one corrected 
for host galaxy extinction. 

The very short duration of optical RS emission in GRB 
06011 IB (about 30 sec in the restframe) raises the question 
of the true fraction of GRBs displaying this type of emission. 
GRB 06011 IB is among the few cases of a burst with bright 
early optical counterpart observed with good temporal resolu- 
tion as well as a 'canonical' Swift/XRT light curve in X-rays 
and provides an excellent test case for the reverse-shock emis- 
sion mechanism in both the energy ranges. 
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